■ii v ui/ 


' % , r : i 

4 A I 

&t>;S L 


PP 

psiM 


jHn ■ 


tiom 


s § r 


D & 


asymptotic expansion 


°F AN IMPULSE FOR AN 


OPTIMAL FINITE BURN 


S. Pines 


(NASA-CR-130243) ASYMPTOTIC EXPANSION OF 
AN IMPULSE FOR AN OPTIMAL FINITE BURN 
(Analytical Mechanics Associates, Inc.) 

11 p H C $3.00 CSCL 22C 



Unclas 

G3/30 17238 ) 


Report No. 72-11 
Contract No. NAS 5-11335 
March 1972 


ANALYTICAL 


mechanics associates 

j0 jer icho turnpike 
Jericho, h. t. 


11753 


Inc. 



SUMMARY 

This report derives closed form expressions for the position and velocity 
of a spacecraft during a finite burn using the method involving the theory of asymp- 
totic expansion of the optimal impulsive solution. The small parameter is given 
by the reciprocal of the mass flow rate. The expansion is given in terms of the 
impulsive solution and holds up through third-order for the velocity and fourth- 
order for the position. 
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DESCRIPTION OF THE OPTIMAL IMPULSIVE SOLUTION 


An optimal-impulse transfer from one trajectory to another satisfies certain 
necessary conditions (Ref. 1). This section will list the pertinent characteristics 
without derivation. It is assumed, in what follows, that the user has available a 
method for achieving the optimal impulsive solution which satisfies these conditions. 

The coasting trajectories, before and after the impulse burn, satisfy the cen- 
tral force field equations 


R = - H 


R 


( 1 ) 


where 


R = f R + gR 
o o 


R = f R + gR 
o o 


(la) 


and f , g , f , and g are scalar functions of the initial conditions and an anomaly 
variable which is related to the time by Kepler's equation (Ref. 2). 
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:r time, wc have 
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The vector A satisfies the variational equations of Eq. (1) 
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The equation for the mass is 


m = - — 
c 


( k , c constant ) 


( 4 ) 


and for the impulse 
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On each of the coasting arcs, the vector A satisfies the transition matrix 
equations of the variational solution of Eq, ( 1 ) 
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( 6 ) 


The necessary optimality condition is given by the characteristic that the 
magnitude of A, X is a maximum at an interval impulse, has a non-positive 
slope at an initial burn, a non-negative slope at a terminal impulse, and that, at 
all impulse times, X is equal to the same maximum value. Furthermore, if X 
exceeds this maximum constant value anywhere on the coasting arcs, the the solu- 
tion is not optimal, and a better solution can be found. 


The Lagrange multipliers, A and A, of the position and velocity state 
must be continuous over the entire arc, except at such times and places where 
constraints are imposed on the state and other irregular instances (Ref. 3). 

We seek a solution of the finite burn solution of the equations 
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R = - 
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in terms of the solution for the impulsive case. 
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EXPANSION OF THE IMPULSE SOLUTION 


Consider the integration of Eq. (7) for a finite time. 


• • p -n r ' k A\ 

R(t) = R(t)-„J -dt + J (--jjdt 

t r t 
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R(t) = R(t o ) + (t-t o )R(t o ) 
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m = m + m( t - 1 ) 
o ' o 


As a zero-order approximation, consider the coast solution of Eq. (1) 
from the same initial conditions 
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Eliminating R(t Q ) and R(t Q ) from Eqs. (8) and (9), we have 
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It is plain from Eq. (10) that, as t-»t , the solution of Eq. (10) approaches 
the impulsive solution given by Eq. (2). We seek a more definitive form of the ex- 
pansion for short finite time. In what follows, let t-t o =T. 


Consider the Taylor series expansion of the difference integrand of the 
gravity term. 
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Since the initial conditions for both coast and burn solutions are the same, the 

2 

first non-vanishing term is the r term. We have 
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where R^ is the transpose of R q and 1 = 3x3 identity matrix. Examination of 

Eq. (lib) shows it to be nonlinear in the vector variable R . Similarly, the next 

3 ' 4 

term, (T ), will be nonlinear in R q and R q and, finally, the T term will be 

nonlinear in the variable A q itself. While we could formally carry out this ex- 
pansion, terms nonlinear in A and A will prove very cumbersome for inversion. 

o o 


In what follows, we restrict the solution to the expansion which is linear in 

the vector variable, A , and its time derivatives. 
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We have the following identity 
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Equation (10) now becomes 
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Before proceeding with the solution, it is timely to point out the time-scale 
difference in the variation of the Lagrange multipliers, A, and the mass flow 
rate. We note that the unit vector in the direction of the thrust rotates relatively 
slowly as compared to the time -variation of the mass. Thus, we take a linear 
variation in thrust direction as a good enough approximation during the burn. We 
have 


_d_ A _ ( Ax A)xA 
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Finally, from the definition of the mass flow rate, we have 
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We define the i^ integral of the reciprocal of the mass 
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In addition, we define the i integral of the first moment of the reciprocal of the 
mass 
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etc. 

The required asymptotic expansion in terms of the reciprocal of the con- 
stant mass flow rate is given by 
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